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Abstract 

We have implemented a Langevin approach for the transport of heavy quarks in the UrQMD hybrid 
model. The UrQMD hybrid approach provides a realistic description of the background medium for the 
evolution of relativistic heavy ion collisions. We have used two different sets of drag and diffusion coefficients, 
' one based on a T-Matrix approach and one based on a resonance model for the elastic scattering of heavy 

9-; quarks within the medium. In case of the resonance model we have investigated the effects of different 

decoupling temperatures of the heavy quarks from the medium, ranging between 130 MeV and 180 McV. 
We present calculations of the nuclear modification factor Raa, as well as of the elliptic flow V2 in Au+Au 
collisions at y/s^N — 200 GeV and Pb+Pb collisions at ^s^n = 2.76 TeV. To make our results comparable 
to experimental data at RHIC and LHC we have implemented a Peterson fragmentation and a quark 
coalescence approach followed by the semileptonic decay of the D- and B-mesons to electrons. We find 
that our results strongly depend on the decoupling temperature and the hadronization mechanism. At a 
decoupling temperature of 130 MeV we reach a good agreement with the measurements at both, RHIC and 
LHC energies, simultaneously for the elliptic flow V2 and the nuclear modification factor Raa- 
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I. INTRODUCTION 



One major goal of ultra- high-energy heavy- ion physics is to recreate the phase of deconfined 
quarks and gluons (the Quark Gluon Plasma, QGP) as it might have existed a few microseconds 
after the Big Bang. Various experimental facilities have been built to explore the properties of this 
QGP experimentally, while on the theory side a multitude of (potential) signatures and properties 
of the QGP have been predicted [3-0] • 

Heavy quarks are an ideal probe for the QGP. They are produced in the beginning of the 
collision in hard processes and therefore probe the created medium during its entire evolution. 
When the system cools down they hadronize, and their decay products can finally be detected. 
By investigating heavy-quark observables we can thus explore the interaction processes within 
the hot and dense medium. Two of the most interesting observables are the nuclear modification 
factor, Raa-, and the elliptic flow, V2- Experimentally, the nuclear modification factor shows a large 
suppression of the open heavy-flavor particles' spectra at high transverse momenta (pr) compared 
to the findings in pp collisions. This indicates a high degree of thermalization also of the heavy 
quarks with the bulk medium consisting of light quarks and gluons and, perhaps at the later stages 
of the fireball evolution, the hot and dense hadron gas. The measured large elliptic flow, V2, of open 
heavy-flavor mesons and the non-photonic single electrons or muons from their semileptonic decay 
underlines this interpretation because it indicates that heavy quarks take part in the collective 
motion of the bulk medium. A quantitative analysis of the degree of thermalization of heavy-quark 
degrees of freedom in terms of the microscopic scattering processes may lead to an understanding of 
the mechanisms underlying the large coupling strength of the QGP and the corresponding transport 
properties. 

In this paper we explore the medium modification of heavy-flavor pr spectra, using a hybrid 
model, consisting of the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model 0, 0| 
and a full (3+l)-dimensional ideal hydrodynamical model jfj|, 0| to simulate the bulk medium. 
The heavy-quark propagation in the medium is described by a relativistic Langevin approach 
0]. Similar studies have recently been performed in a thermal fireball model with a combined 
coalescence- fragmentation approach in an ideal hydrodynamics model with a lattice-QCD 



EoS [mjl6|, in a model from Kolb and Heinz [12|, in the BAMPS model 113,1131, the MARTINI 



model [20| as well as in further studies and model comparisons [2114251] . 

The use of the UrQMD hybrid model provides a major step forward as compared to simplified 
expanding fireball models employed so far. It provides a realistic and well established background, 
including event-by-event fluctuations and has been shown to describe well many collective prop- 
erties of relativistic heavy- ion collisions. For the heavy-quark propagation we apply a Langevin 
approach. Within this framework we investigate the effects of using different drag and diffusion co- 
efficients and different freeze-out temperatures of heavy flavors on the heavy-quark observables and 
compare the results with the experimental data from the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC). 



II. THE URQMD HYBRID MODEL 

To extract information on the interaction of heavy quarks with the medium one ideally applies 
a well tested model for the (collective) dynamics of the bulk matter. In heavy-ion collisions the 
medium is by no means homogeneous. Rather it is a locally and event-by-event fluctuating, fast 
expanding system. In our calculation we employ the state of the art UrQMD hybrid model for 
the description of the expanding background. This model has been developed in the past years to 
combine the advantages of hadronic transport theory and ideal fluid dynamics J26[. To account 
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for the non-equilibrium dynamics in the very early stage of the collision in the hybrid model, the 
UrQMD cascade 0, @] is used to calculate the initial states of the heavy ion collisions each to 



be used in a subsequent hydrodynamical evolution 271 ] . The transition from the UrQMD initial 
state and the hydrodynamical evolution takes place at a time t star t = 2i?/ \J 7^ M — 1 i.e., after the 
two Lorentz-contracted nuclei have passed through each other (7cm is the center-of-mass-frame 
Lorentz factor, and R is the radius of the nucleus). The energy, baryon number, and momenta 
of all particles within UrQMD are mapped onto a spatial grid for the hydrodynamic evolution 
including event-by-event fluctuations. The full (3+l)-dimensional ideal hydrodynamic evolution 
is performed using the SHASTA algorithm 0,0). We solve the equations for the conservation of 
energy and momentum and for the conservation of the baryonic charge. With T^ v denoting the 
relativistic energy-momentum tensor the corresponding equations read 



= 0, (1) 



and for the baryon four-current 



= 0. (2) 
To transfer all particles back into the UrQMD model, an approximate iso-eigentime transition is 



chosen (see 28| for details). Here, we apply the Cooper- Frye prescription [29( and transform to 



particle degrees of freedom via 



dN 



Here da^ = (d 3 x, 0,0,0)) is the hypersurface normal. In Eq. ([3]) f(x,p) are the Bose- and Fermi- 
distribution functions and gi the degeneracy factors for the different particle species. After the 
"particlization" the evolution proceeds in the hadronic cascade (UrQMD), where final re-scatterings 
and decays are calculated until all interactions cease and the system decouples. 

A more detailed description of the hybrid model including parameter tests and results can be 



found in 26]. A comparison to the results employing the non-approximated hypersurface can be 
found in 0. 



III. HEAVY-QUARK DIFFUSION 

The diffusion of a "heavy particle" in a medium consisting of "light particles" can be described 
with a Fokker-Planck equation 0, I21I . 31-35)]. Here one approximates the collision term of the 



corresponding Boltzmann equation, which in turn can be mapped into an equivalent stochastic 
Langevin equation. 

A. Relativistic Langevin approach 

In the relativistic realm such a Langevin process reads 

dxj = -^dt, 

E r- (4) 

dpj = —Tpjdt + VdtCjkPk- 

Here dt is the time step in the Langevin calculation, dxj and dpj are the coordinate and momentum 
changes in each time-step, E = y m 2 + p 2 , and T is the drag or friction coefficient. The covariance 
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matrix, Cjk, of the fluctuating force is related to the diffusion coefficients, as we shall see below. 
Both r and Cjk dependent on (t,x,p) and are defined in the (local) rest-frame of the fluid. The 
Pk are Gaussian- normal distributed random variables. Their distribution function reads 

with p = (pi, p2, /O3). The fluctuating force F- thus obeys 

(if } (i)) = 0, (Ff(i)ifV)) = C fl C u 5{t - t'). (6) 

It is important to note that with these specifications the random process is not yet uniquely 
determined since one has to specify, at which momentum argument the covariance matrix Cjk has 
to be taken to define the stochastic time integral in Q. Thus, we set 

C jk = Cjk{t,x,p + t&p). (7) 

For £ = 0, ^ = 1/2, and £ = 1 the corresponding Langevin processes are called the pre-point 
Ito, the mid-point Stratonovic-Fisk, and the post-point Ito (or Hanggi-Klimontovich) realization, 
respectively [361 ], 

According to ([!]) and ©, for a given value of £ in (|7|) the average of an arbitrary observable 
g{x,p) obeys the time-evolution 

w . + d,, P + d P) - g( *, P )> (;;«■*■:;« (-r n + 

Here all momentum arguments of the drag and diffusion coefficients have to be taken at p. From ([8]) 
it follows immediately that the time evolution of the phase-space distribution function fQ(t,x,p) 
of heavy quarks is given by the Fokker-Planck equation, 



dfg Pj df Q = _d_ 
dt E dxj dpj 



Thus, the usual drag and diffusion coefficients for an isotropic medium are related to the pertinent 
parameters in the Langevin process by 

BC 

Apj=T P j-iC lk -^, (10) 
Cj k = v^^i + VlBlPfk, (11) 
with 4 = ^r> i% = *i*- P -f- (12) 

case of a homogeneous static background ("heat bath"), the stationary limit should become 
a Boltzmann- Jiittner distribution with the temperature of the "heat bath" . Thus, one typically 
adjusts the drag coefficient by choosing the longitudinal diffusion coefficient, Bi, in (fTT|) such as to 



1 In numerical studies it has turned out that drag and diffusion coefficients as obtained from microscopic models 
usually do not lead to the expected long-time stationary limit of the phase-space distribution for the heavy particles 
when diffusing in an equilibrated background medium. 
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satisfy this asymptotic equilibration condition, leading to dissipation-fluctuation relations between 
this diffusion coefficient and the drag coefficient @, 21]. 

It turns out that for Bq = B\ = D(E) and a homogeneous background medium the Boltzmann- 
Jiittner distribution, 



/q cci) (p) = exp y-ipj » with E = Vp 2 + m2 , 



(13) 



becomes a solution of the corresponding stationary Fokker-Planck equation, if the dissipation- 
fluctuation relation 



T(E)ET-D(E)+T(1-£)D'(E) = 0, (14) 

is fulfilled. A straightforward way to achieve the correct asymptotic equilibrium distribution within 
a relativistic Langevin simulation is to set £ = 1 (i.e., using the post-point Ito realization). This 
reduces (fl4l) to 

D{E) = T(E)ET. (15) 

For applications to heavy-ion collisions we use T and Bq from underlying microscopic models for 
heavy-quark scattering with light quarks and gluons as detailed below and adjust the longitudinal 
diffusion coefficient to 



B x = VET. (16) 

So far we have defined our Langevin process with respect to the (local) rest frame of the background 
medium. For a medium with collective flow, one has to evaluate the time-step in the local rest- 
frame and boost-back to the computational frame. For a closer look on the post-point description 
see section IVH Al 

For the heavy-quark propagation in the Langevin model we also need transport coefficients. In 
this work these drag and diffusion coefficients are obtained from two non-perturbative models for 
elastic heavy-quark scattering, a resonance model, where the existence of D-mesons and B-mesons 
in the QGP phase is assumed, as well as a T-Matrix approach in which quark-antiquark potentials 
are used for the calculation of the coefficients in the QGP. They are described in detail in Sec. 
IVII Bl and are shown in Fig. Q]as function of the three- momentum \p\ at T = 180 MeV and in Fig. 
[2] as function of the temperature at a fixed three-momentum of \p\ = 0. 



B. Implementation of the Langevin simulation into the UrQMD-hybrid model 



For the present study, charm production and propagation is evaluated perturbatively on the 
time-dependent background generated by UrQMD/Hybrid. To model a fluctuating and space-time 
dependent Glauber-initial state geometry, we perform a first UrQMD run with elastic 0° scatterings 
between the colliding nuclei and save the nucleon-nucleon collision space-time coordinates. These 
coordinates are used in a second, full UrQMD run as (possible) production coordinates for the 
heavy quarks. 

As momentum distribution for the initially produced charm quarks at ^/snn = 200 GeV we 

use 

_JL_ = J^ML (17) 



2TTp T dp T (l + A 2 -p 2 T ) 



5 



0.18 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




T=180 MeV 



charm resonance model 
charm T-Matrix approach 
bottom resonance model 
bottom T-Matrix approach 





0.18 




0.16 




0.14 




0.12 


E 




CVl 


0.1 


> 




CD 




CD 


0.08 


o 
CD 






0.06 




0.04 




0.02 








T=180 MeV 


charm resonance model 


charm T-Matrix approach . 




bottom resonance model 




bottom T-Matrix approach 







p[GeV] 



3 4 
p[GeV] 



FIG. 1. (Color online) Drag (left) and diffusion (right) coefficients in the resonance model and the T- 
Matrix approach for charm and bottom quarks. The plot shows the dependence of the coefficients on the 
three-momentum \p\ at a fixed temperature of T = 180 MeV. 



10 



p=0 GeV 



E 
< 



0.1 



charm resonance model 
charm T-Matrix approach 
bottom resonance model 
bottom T-Matrix approach 



0.01 



100 



200 300 
T[MeV] 



400 



500 



> 

CD 

a 

o 
CD 



10 



0.1 r 



0.01 



p=0 GeV 


charm resonance model 

charm T-Matrix approach 




bottom resonance model 




bottom T-Matrix approach 















100 



200 300 
T[MeV] 



400 



500 



FIG. 2. (Color online) Drag (left) and diffusion (right) coefficients in the resonance model and the T- 
Matrix approach for charm and bottom quarks. The plot shows the dependence of the coefficients on the 
temperature at a fixed three-momentum |p| = 0. The T-Matrix coefficients are calculated between f 80 MeV 
and 360 MeV only. 



with Ai = 0.5, A2 = 0.1471 and A3 = 21 and for bottom quarks 



2np T dpT [A 1+ pl) M ' 

with A\ = 57.74 and A2 = 5.04. These distributions are taken from 0, [37j and are shown in Fig. 
[3l They are obtained by using tuned c-quark spectra from PYTHIA. Their pertinent semileptonic 
single-electron decay spectra account for pp and dAu measurements by STAR up to 4 GeV. The 
missing part at higher pr is then supplemented by B-meson contributions. 

Starting with these charm- and bottom-quark distributions as initial conditions we perform, as 
soon as the hydrodynamics start condition is fulfilled, an Ito post-point time-step of our Langevin 
simulation as described in Sec. IIII Al at each time-step of the hydro dynamical evolution. 

We use the cell velocities, cell temperatures, the length of the time-step and the 7-factor of the 
cells to calculate the momentum transfer, propagating all heavy quarks independently. For the 
Langevin transport we use the drag and diffusion coefficients obtained from the resonance model 
or T-Matrix approach as described in Sec. IVIIBI 

To analyze the sensitivity of Raa and especially V2 on the decoupling time of the heavy fla- 
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FIG. 3. (Color online) Fits of D- and D* -meson px spectra in 200 AGeV d-Au collisions at RHIC with a 
modified PYTHIA simulation (left panel) and the corresponding non-photonic single-electron pr spectra in 
p-p and d- Au collisions (taken from j38| ) . The missing yield of high-pT electrons is fitted with the analogous 
B-meson decay spectra, thus fixing the bottom-charm ratio at a b i/cr C c — 4.9 • 10~ 3 . 



vors from the medium we vary the decoupling temperatures between 130 MeV and 180 MeV (for 
the resonance model) and extrapolate the corresponding transport coefficients smoothly into the 
hadronic phase. This assumption of a smooth transition of the transport coefficients in the tran- 
sition from the partonic description above and the hadronic one below T c has been verified, using 
an effective model for open- heavy-flavor interactions in a hadronic medium in 15|, |39( . 

Our approach provides us with the heavy-quark momentum distribution. We include a 
hadronization mechanism for open-heavy-flavor mesons (D and B mesons). Since non-photonic 
single electrons are usually measured in experiments, we perform a semileptonic decay into elec- 
trons as final step to compare to data. In addition we also provide D- and B-meson results for 
direct comparisons to the upcoming direct D /B measurements by the STAR Heavy Flavor Tracker 
(HFT). These results are shown in Sec. IVIICl 



IV. RESULTS AT RHIC ENERGIES 

A. Elliptic flow v-2 and nuclear modification factor Raa with fragmentation 

Fig. H] presents the elliptic flow, V2, of charm and bottom quarks from Au+Au collisions at 
y/ s NN = 200 GeV in the centrality range a jo tot =20%-40% applying a rapidity cut of \y\ < 0.35. 

For our calculation using the drag and diffusion coefficients of the T-Matrix model we use a 
decoupling temperature of 180 MeV, while with the resonance model we show results for decoupling 
temperatures of 130 MeV, 150 MeV and 180 MeV. 

As one can clearly see, the elliptic flow, V2, of bottom quarks (dashed lines) is much smaller 
compared to that of the charm quarks (solid lines) due to their larger mass. Furthermore the 
use of the coefficients from the T-Matrix model compared with those from the resonance model 
shows that both calculations are in reasonable agreement. The elliptic flow of the charm quarks 
is nevertheless somewhat lower for the T-Matrix model than for the resonance model. When 
decreasing the decoupling temperature the flow clearly increases. Thus, we conclude that the late 
phase of the heavy-ion collision may have considerable influence on the heavy-flavor elliptic flow 
although the drag and diffusion coefficients become small in the late stages of the fireball evolution. 

Moreover the, V2, is shifted towards higher px for lower decoupling temperatures. This effect 
is due to the increased radial velocity of the medium, which is in case of an developed elliptic flow 
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larger in x than in y direction. Consequently there is a depletion of particles with high v x in the 
low pt region and smaller elliptic flow. This effect is more important for heavier particles and a 
larger radial flow 40j, [41 ] . 

To compare our calculations with data on non-photonic electrons from RHIC we perform a 
Peterson fragmentation of the charm and bottom quarks to D-mesons and B-mesons using the 
fragmentation function from 421 ]. 



N 



z[1-(1/z)-€q/{1-z)]' 



where N is a normalization constant, z the relative-momentum fraction obtained in the fragmen- 
tation of the heavy quark and eg = 0.05(0.005) for charm (bottom) quarks. After hadronization 



we use PYTHIA routines for the semileptonic decay to electrons [42 

Fig. [5] shows our results for the V2 for single-electrons in comparison to data from the PHENIX 
collaboration. 
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FIG. 4. (Color online) Left: Elliptic flow, t>2, of heavy quarks in Au+Au collisions at ^/JWn = 200 GeV. 
We use a rapidity cut of \y\ < 0.35. The solid lines depict the charm quarks while the dashed lines depict 
the bottom quarks. 

Right: Elliptic flow, w 2 , of electrons from heavy- meson decays using Peterson fragmentation to D/B mesons 
and subsequent decay into electrons in Au+Au collisions at ^/snn = 200 GeV. We use a rapidity cut of 
\y\ < 0.35. Data are from jlij ]. 

Again we clearly observe the importance of the late phase of the collision. The depletion effect 
at low pt described before is clearly visible. The decrease of the elliptic flow at high p? is due to 
the increasing fraction of electrons from bottom decays, which have a lower v% as seen in Fig. [H 
The calculated flow in the setup with the Peterson fragmentation is too small compared to the 
PHENIX data. 

The corresponding nuclear modification factor, Raa, for heavy quarks is shown in Fig. [5j Again 
we present results for Au+Au collisions at ^Jsnn = 200 GeV in the centrality range 20%-40%. 
The quenching for charm quarks is, as expected, much stronger than for bottom quark^l. While 
for bottom quarks the suppression at high px is moderate, Raa may drop to 20-30% for charm 
quarks. The influence of the medium is, as already seen in our flow calculations, larger for a lower 
decoupling temperature underlining the importance of the late phase of the collision. Fig. [5] shows 
the comparison of our non-photonic-electron Raa to the data taken by the PHENIX collaboration. 



2 However, recent preliminary PHENIX data presented at QM 2012 seem to suggest the contrary. 
(Talk by T. Sakaguchi at QM 2012, data not published yet) 
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FIG. 5. (Color online) Left: Raa of heavy quarks in Au+Au collisions at ^Jsnn = 200 GeV. We use a 
rapidity cut of \y\ < 0.35. The solid lines depict the charm quarks while the dashed lines depict the bottom 
quarks. 

Right: Raa of electrons from heavy quark decays in Au+Au collisions at y/spjN = 200 GeV compared to 
RHIC data [45J]. We use a rapidity cut of \y\ < 0.35. The high-p^ suppression turns out to be too strong 
compared with the data. 



The nuclear modification factor drops quite rapidly and stabilizes at about px ~> 2 GeV. Around 
Pt ~ 2 GeV it is significantly below the PHENIX data. For higher px the calculated Raa ap- 
proaches the measured data, especially for low decoupling temperatures. This effect is due to the 
increasing flow of the heavy-flavor particles with decreasing decoupling temperature, which pushes 
\ow-pt heavy-flavor particles towards higher pt bins. 



B. Elliptic flow V2 and nuclear modification factor Raa using a k factor 

In the previous section we learned that the elliptic flow of the heavy quarks in the calculation 
with fragmentation is too small compared to experimental data. One possibility to improve on 
this problem may be to multiply the drag and diffusion coefficients with a "fc factor". Therefore 
we have performed the same calculations as in the last section but using a k factor of 3. 

As we see in Fig. [6] the elliptic flow increases considerably due to the stronger coupling of the 
heavy quarks to the hot medium. The results after performing the Peterson fragmentation and 
the subsequent decays to electrons are shown in Fig. [6J The elliptic flow is now comparable to 
the data, especially when using a low decoupling temperature of 130 MeV. Only at low pt we 
underestimate the flow due to the depletion effect described above. 

Our results for the nuclear modification factor, Raa-, are depicted in Fig. [3 The quenching is 
much stronger than for the calculation without a k factor. Fig. [7] shows the results for electrons. 
The suppression of non-photonic electrons at high px is also stronger than for the calculation 
without a k factor. 

We conclude that the use of a k factor can improve the description of the elliptic flow. However, 
it is not possible to reach a consistent simultaneous description of both, the elliptic flow and the 
nuclear modification factor using the same k factor. 



C. Elliptic flow «2 and nuclear modification factor Raa using Coalescence 

Instead of describing heavy quark hadronization by Peterson fragmentation (and / or and addi- 
tional k-factor, as discussed above) one can alternatively apply a quark coalescence approach for 
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FIG. 6. (Color online) Left: Elliptic flow, V2, of heavy quarks in Au+Au collisions at -y/iivjv = 200 GeV 
employing a k factor of 3. We use a rapidity cut of \y\ < 0.35. The solid lines depict the charm quarks while 
the dashed lines depict the bottom quarks. 

Right: Elliptic flow, V2, of electrons from heavy quark decays in Au+Au collisions at ^/sWn = 200 GeV 
employing a k factor of 3. We use a rapidity cut of \y\ < 0.35. The flow in our calculation using a k factor 
is comparable to data 45 1. 
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FIG. 7. (Color online) Left: Raa of electrons from heavy-flavor decays in Au+Au collisions at ^snn = 
200 GeV employing a k factor of 3. We use a rapidity cut of \y\ < 0.35. The solid lines show the results for 
charm quarks and the dashed ones for bottom quarks. 

Right: V2 of electrons from heavy-flavor decays in Au+Au collisions at y^s^T/v = 200 GeV employing a k 
factor of 3. We use a rapidity cut of \y\ < 0.35. As expected the medium modification is stronger than 
without a k factor. Data are taken from 1451 . 



D- and B-meson production. To implement this coalescence we perform the Langevin calculation 
until the decoupling temperature is reached. Subsequently we coalesce a light quark with a heavy 
quark. As the light quarks constitute the medium propagated by hydrodynamics, the average 
velocities of the light quarks can be (on average) approximated by the flow-velocities of the hydro 
cells. The mass of the light quarks is assumed to be 370 MeV so that the D-meson mass becomes 
1.87 GeV when the masses of the light quarks and the charm quarks (1.5 GeV) are added. Since 
we assume the light quarks to have the same mass when coalescing with bottom quarks (4.5 GeV), 
the B-mesons obtain a mass of 4.87 GeV. 

The differences of the flow and the spectra of D- and B-mesons when comparing Peterson 
fragmentation (without fc-factor) to the coalescence model is shown in Fig. [8J These calculations 
are performed employing a decoupling temperature of 150 MeV. 

As compared to the fragmentation case, the elliptic flow reaches higher values at high px due 
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FIG. 8. (Color online) Elliptic flow (left) and Raa (right) of D mesons (solid lines) and B mesons (dashed 
lines) in Au+Au collisions at y/s^N = 200 GeV. We use a rapidity cut of \y\ < 0.35. A comparison of a 
Peterson fragmentation and a coalescence with light quarks is shown. For the drag and diffusion coefficients 
we use the resonance model with a decoupling temperature of 150 MeV. 



to the coalescence. Also the depletion effect described before is more pronounced. Regarding the 
nuclear modification factor, Fig. [SJ the difference of Peterson fragmentation and the coalescence 
model is even larger. The push of low-pr particles to higher pr is stronger in case of the coa- 
lescence model, while the suppression of heavy mesons at high pj- is stronger in case of Peterson 
fragmentation. 

Again we perform a decay to electrons using PYTHIA to compare to experimental measurements 
from the PHENIX collaboration. Fig. [9] (left) shows our results for V2- Due to the coalescence 
the elliptic flow is strongly increased compared to the previous calculation using the Peterson 
fragmentation. This higher flow is due to the momentum kick of the light quarks in the recombi- 
nation process, which provides additional flow from the medium. For a decoupling temperature of 
130 MeV we obtain a reasonable agreement with the experimental data. 

In Fig. [9] (right) the nuclear modification factor for non-photonic single electrons is depicted. 
Also here we obtain a good agreement with the data. Especially at moderate pt ~ 2 GeV the 
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FIG. 9. (Color online) Elliptic flow V2 (left) and nuclear modification factor Raa (right) of electrons from 
heavy quark decays in Au+Au collisions at ^/snn = 200 GeV using a coalescence mechanism. We use a 
rapidity cut of |y| < 0.35. For a decoupling temperature of 130 MeV we get a reasonable agreement to data 

calculation has strongly improved. The coalescence mechanism pushes the heavy quarks to higher 
Pt- As seen before we obtain the best agreement to data for rather low decoupling temperatures. 
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In conclusion we observe that the coalescence mechanism is required to describe experimental 
data with our Langevin model. Only with the coalescence model one is able to describe both Raa 
and V2 consistently in the present model. 



D. Dependence of the medium modification on the equation of state 

The heavy-flavor-flow observables in Langevin simulations are quite sensitive to the used descrip- 



tion of the background medium |24l ]. To examine this issue somewhat further, we have performed 
our calculations also using different equations of state that are implemented in the UrQMD hybrid 
model. Our results for different equations of state for the drag and diffusion coefficients of the 
resonance model with a decoupling temperature of 150 MeV are shown in Fig. [10] for the elliptic 
flow V2 and for the nuclear modification factor Raa- 

The equation of state we have been using for all results in the previous sections is the chiral EoS. 
It is constructed by matching a state of the art hadronic chiral model to a mean field description 
of the deconfined phase. The deconfinement transition in this approach is included by the means 
of an effective Polyakov Loop potential, coupling to the free quarks. It has been shown in [iH ] that 
the chiral EoS gives a reasonable description of lattice QCD thermodynamics at fis = and can 
be extended to finite baryon densities. The Hadron resonance gas EoS resembles the active degrees 
of freedom that are also included in the UrQMD transport approach, namely most hadronic states 
and their resonances. The Bag model EoS |?J follows from matching a Walecka type hadronic 
model to massless quarks and gluons via a Maxwell construction. It exhibits a strong first order 
phase transition for all values of hb- 




FIG. 10. (Color online) Elliptic flow t>2 (left) and nuclear modification factor Raa (right) of heavy quarks 
in Au+Au collisions at ^/sjvat = 200 GeV. We use a rapidity cut of \y\ < 0.35. Different equations of state 
are compared. 

As one sees clearly the influence on the mediums evolution as seen through the heavy quarks 
for this set of equations of state is very small. 



V. RESULTS AT LHC ENERGIES 

In the previous sections we found that we reach the best agreement to experimental PHENIX 
data when using the Resonance model applying a decoupling temperature of 130 MeV and using 
quark coalescence as hadronization mechanism. Now we apply the same description also at LHC 
energies (y / s/v r /v = 2.76 TeV). The momentum distribution for the initially produced charm quarks 
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at LHC is obtained from a fit to PYTHIA calculations. The fit function we use is 



d*pr {i + A l -{pl) M )M 

with the coefficients A\ = 0.136, A 2 = 2.055 and A 3 = 2.862. 

We have performed our calculations in Pb+Pb collisions at ^snn = 2.76 TeV in a centrality 
range of 30%-50%. The analysis is done in a rapidity cut of \y\ < 0.35 in line with the ALICE 
data. 
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FIG. 11. (Color online) Left: Flow v-i of D-mesons in Pb+Pb collisions at *J~sWn = 2.76 TeV compared to 
data from the ALICE experiment. (Talk by Z. Conesa del Valle at QM 2012, data not published yet.) A 
rapidity cut of \y\ < 0.35 is employed. 

Right: Raa of D-mesons in Pb Pb collisions at v / s]v7V = 2.76 TeV compared to experimental data from 



ALICE [47J. A rapidity cut of \y\ < 0.35 is employed. 



Fig- LTD (left) depicts our results for the elliptic flow compared to ALICE measurements. The D- 
meson v 2 exhibits a strong increase and reaches a maximum at about px = 3 GeV with v 2 ~ 19%. 
The agreement between the ALICE measurements of D°- and D + -mesons and our calculation is 
quite satisfactory. 

A complementary view on the drag and diffusion coefficients is provided by the nuclear sup- 
pression factor Raa- Figure [TTJ (right) shows the calculated nuclear modification factor Raa of 
D-mesons at LHC. In line with the experimental data the simulation is done for a more central bin 
of cr/otot = 0%-20%. We find a maximum of the Raa at about pr = 2 GeV followed by a sharp 
decline to an Raa of about 0.2 at high pr. As we can see we can describe the data at medium pt 
well but over-predict them at low pt bins. 



VI. SUMMARY 



In this paper we have investigated the medium modification of heavy-quark pt spectra in the 
hot medium created in heavy-ion collisions at RHIC and LHC energies on the basis of a Langevin 
simulation coupled to the UrQMD hybrid model. The aim of this study was to find a consistent 
description for both the elliptic flow, v 2 , and the nuclear modification factor, Raa, with a realistic 
dynamical description of the background medium. We have used two different sets of drag and 
diffusion coefficients, based on a T-Matrix approach and a resonance-scattering model for the 
elastic scattering of heavy quarks with light quarks and antiquarks. Both sets of coefficients lead 
to similar results for the heavy-flavor observables. 
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In the first part of our analysis we have used a fragmentation mechanism, the Peterson frag- 
mentation, to describe the hadronization of heavy quarks to open-heavy-flavor mesons. We have 
found a low elliptic flow and a too strong heavy- flavor suppression at high pr- Subsequently we 
have explored how a k factor for the drag and diffusion coefficients would influence the results. We 
found that with k = 3, the description of vi is improved, but has lead to an even larger suppression 
of the nuclear modification factor Raa, as expected. We conclude that a combination of fragmen- 
tation and a Langevin simulation with a fc-factor in the transport coefficient does not allow for 
a consistent description of the data on non-photonic single electron spectra in Au+Au collisions 
(y/smr = 200 GeV) at RHIC. 

To overcome this problem we have used a coalescence approach to heavy-quark hadronization 
to open-heavy-flavor mesons instead of the fragmentation. The coalescence mechanism allows for a 
consistent description of both vi and Raa- We have performed the simulations, assuming different 
decoupling temperatures of the heavy quarks from the medium, and found that the late phase of 
the collision can have a considerable effect on the heavy-quark observables. Within our study we 
find the best agreement to experimental data using a low decoupling temperature of 130 MeV. In 
Sec. IIVDI we have also addressed the sensitivity of the heavy-flavor observables to the assumed 
equation of state of the strongly interacting medium. Here we find that our results are insensitive 
to variations of the particular equation of state used in UrQMD's hydrodynamic model. 

Finally we also explored the medium modification in our model at LHC energies. Here we could 
reach a good agreement to data for the elliptic flow V2 of D-mesons. For the nuclear modification 
factor Raa we reach a good agreement at medium pr, but seem to miss the data at low pt bins. 

First measurements with the STAR Heavy Flavor Tracker (HFT) are scheduled for the year 2014. 
It will provide new, complementary measurements in heavy ion collisions at ^/snn = 200 GeV. 
The HFT will enable direct identification of heavy flavor meson decays like D° — > K~ir + or 
— > K~ir + K + . This is supposed to lead to better vi measurements down to very low pt and a 
better understanding of the energy loss of heavy quarks in the medium. Especially it will provide 
us with identified D-meson spectra which will enable us to compare our heavy-meson results to 
data separately for D- and B-mesons and therefore to get further insights on the hadronization 
mechanism. 
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VII. APPENDIX 



A. Post-point Ito realization 

Since the phase-space distribution of relativistic particles is a scalar [4§], the proper equilibrium 
limit is given by the corresponding boosted Boltzmann-Jiittner phase-space distribution, 



/^ocex P (-^j, (20) 

where u(t, x) is the four-velocity field of the medium and p the (on-shell) four-momentum of 
the heavy quark in the local rest-frame. It can be shown analytically, and we have numerically 
checked, that for obeying this constraint, one has to apply the post-point prescription, £ = 1, 
strictly only to the momentum argument of the covariance matrix, Cj^ as given in ([7]) and not to 
the corresponding coefficients originating from the Lorentz transformation of the time step dt with 
respect to the laboratory frame (bare symbols) to the one in the local rest-frame of the heat bath 
(starred symbols), i.e., in the transformation prescription for the time interval, 

, m m E E , . 

dt* = — dr = — -dt = dt, 21 

E* E* m p ■ u 

one has to use the heavy-quark momenta at time t without a post-point update rule. Here, dr 
denotes the scalar "proper-time" interval of the heavy quark, corresponding to the given time 



interval, dt, with respect to the laboratory frame [241 ]. 



B. Drag and diffusion coefficients 

We use two non-perturbative models for elastic heavy-quark scattering in the quark-gluon 
plasma to evaluate the drag and diffusion coefficients for the Langevin simulation of heavy-quark 
diffusion. 

The resonance model is based on heavy-quark effective theory (HQET) and chiral symmetry 
in the light-quark sector (33|. Motivated by the finding in lattice-QCD calculations that hadron- 
like bound states and/or resonances might survive the phase transition in both the light-quark 
sector (e.g., p mesons) and heavy quarkonia (e.g., J/ip), in this model we assume the existence of 
open- heavy- heavy- flavor meson resonances like the D and B mesons. 

In the T-Matrix approach static in-medium quark-antiquark potentials from lattice QCD are 
used as scattering kernels in a Bruckner like T-matrix approach to calculate the scattering-matrix 
elements for elastic scattering of heavy quarks with light quarks and antiquarks 0] . 

The heavy-light quark resonance modelf^ is based on the Lagrangian, 



& Dcq =K + ^c° q ~ im^-c - qi^Uf-c + h.c 
- G v (q^$* l±i c - g 7 5 7 ^ lAt l±i c + h.c. 



(22) 



and an equivalent one for bottom quarks. Here v denotes the heavy-quark four- velocity. The free 
part of the Lagrangian is given by 

_^ c ° 9 = c(i0 - m c )c + qi0q, 

= (d^)(d^) + (fyVtya^S) - m|(*t$ + $*t $ * } (23) 
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where <E> and $q are pseudo-scalar and scalar meson fields (corresponding to D and Dq mesons). 
Based on chiral symmetry, restored in the QGP phase, we also assume the existence of mass 
degenerate chiral-partner states. Further from heavy-quark effective symmetry one expects spin 
independence for both the masses, mj = my, and the coupling constants, Gs = Gy. For the 
strange-quark states we take into account only the pseudo-scalar and vector states (D s and D*, 
respectively). 

The D-meson propagators are dressed with the corresponding one-loop self energy. Assuming 
charm- and bottom-quark masses of m c = 1.5 GeV and mj, = 4.5 GeV, we adjust the masses of the 
physical D-meson- like resonances to mo = 2 GeV and tub = 5 GeV, in approximate agreement 



with the T-matrix models of heavy- light quark interactions in 49|, |50(] . The coupling constant is 
chosen such as to obtain resonance widths of Td,b = 0.75 GeV. 

With these propagators the elastic Qq- and Qg-scattering matrix elements are calculated and 
used for evaluation of the pertinent drag and diffusion coefficients for the heavy quarks, using (|29|) 
and (|30p . It turns out that particularly the s-channel processes through a D /B-meson like resonance 
provide a large efficiency for heavy-quark diffusion compared to the pQCD cross sections for the 
same elastic scattering processes. This results in charm-quark equilibration times tZ, = 2-10 fm/c. 

In order to justify the formation of D- and B-meson like resonances above T c , in [9j a Brueckner- 
like in-medium T-matrix approach has been used for the description of elastic heavy-light-quark 
scattering in the QGP. After a three-dimensional reduction to a Lippmann-Schwinger equation, 
including a Breit correction, in-medium heavy-quark potentials from 1QCD have been employed as 
the scattering kernels. As an upper limit of the interaction strength within such an approach, the 
internal-energy potential, 

U(r,T) = F(r,T)-T^f±, (24) 

has been used, where F is the free-energy potential from the lattice calculation. We take into 
account also the complete set of Qq color states, assuming Casimir scaling of the corresponding 
potentials, 

V 8 = ~V 1 , Kf=ivi, V 6 = -\v 1 . (25) 
After a partial-wave decomposition the Lippmann-Schwinger equation, 

T a j(E;q',q) = V a ,i(q',q) 

+ ljdk k 2 V a ^(q', k)G Qq (E, k) (26) 

xT n , ( (fi;M)[l-/ f K°)-/fK)], 

for the partial-wave components of each color channel, a, has been solved for the S- and P-wave 
components. Here, E is the center-of-momentum (CM) energy of the heavy-light quark system, q 
and q' the momenta of the heavy and light quark, and 

G qQ (E, k) = j= (27) 

qQ £-K + i£?)-K« + iX?) 

the corresponding two-particle propagator in the CM frame. It has been checked that the quasi- 
particle widths of Tf Q = 2S|' Q = 200 MeV are consistent with a previous similar Br"uckner 



calculation [5l|] for the light quarks and with the heavy-quark self-energies with the T-matrix 
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solution of (|26p . The relation with the invariant scattering-matrix elements in (|29|) is then given 
by 



^2\M\ 2 = i^(s - m\ + m|) 2 (s - - m 2 ) 2 



(28) 

=0 (s)| 2 + 3|T aifc= i(s)cos6> cm | 2 ) . 

a 

The relation of elastic heavy-quark-scattering matrix elements with the drag and diffusion coeffi- 
cients in the Langevin approach is given by integrals of the form 

, , m = 1 f d 3 q f d 3 p> f dV 1 v^|.,|2 

^ [P )} ~2u p 7 R 3 2£(g) (2vr)3 7 M 3 2£(p') (2vr)3 7 M 3 2E(q>) (2tt)3 7Q ' ' ( 29 ) 

x(27r)V 4 )(p + (? -p / -g')/ g , 9 (q)^(p') • 

Here, the integrations run over the three momenta of the incoming light quark or gluon and the 
momenta of the outgoing particles. The sum over the matrix element is taken over the spin and 
color degrees of freedom of both the incoming and outgoing particles; 7q = 6 is the corresponding 
spin-color degeneracy factor for the incoming heavy quark, and f qg stands for the Boltzmann 
distribution function for the incoming light quark or gluon. In this notation, the drag and diffusion 
coefficients are given by 

B (P) = \ (p' 2 - ^f- ) , (30) 



5i(p) = ^(4^-Vp + P 



1 /(PPL 
2\ p2 

For both approaches we also include the leading-order perturbative QCD cross sections for 
elastic gluon heavy-quark scattering 52| , including a Debye screening mass rriDg = gT in the gluon 
propagators, thus taming the i-channel singularities in the matrix elements. The strong-coupling 
constant is chosen as a s = g 2 /(4tt) = 0.4. 



17 



C. Underlying D- and B-meson spectra before semi-leptonic decays 



The heavy flavor electron spectra at RHIC originate from D- and B-meson decays. These D- 
and B-meson spectra are obtained from our heavy quark calculations applying a fragmentation or 
a coalescence mechanism. They are displayed in Fig. [12] for the case of the Peterson fragmentation 
without using a /c-factor, in Fig. [13] for the case of the Peterson fragmentation applying a /c-factor 
of 3 and finally for the case of using a coalescence mechanism (Fig. [bl|h 

These spectra can act as a prediction for future D- and B-meson measurements at RHIC energies. 
On the one hand they allow for a comparison of our hadronization mechanisms to experimantal 
data and on the other hand for a comparison of the decay to heavy flavor electrons performed using 
PYTHIA. 




FIG. 12. (Color online) Elliptic flow V2 (left) and nuclear modification factor Raa (right) of D- and B- 
mesons using Peterson fragmentation in Au+Au collisions at ^/s^n = 200 GeV. We use a rapidity cut of 
\y\ < 0.35. 




FIG. 13. (Color online) Elliptic flow V2 (left) and nuclear modification factor Raa (right) of D- and B- 
mesons using Peterson fragmentation and a k-factor of 3 in Au+Au collisions at ^/sjvat = 200 GcV. We use 
a rapidity cut of \y\ < 0.35. 
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FIG. 14. (Color online) Elliptic flow V2 (left) and nuclear modification factor Raa (right) of D- and B-mesons 
using coalescence in Au+Au collisions at ^/snn = 200 GeV. We use a rapidity cut of \y\ < 0.35. 
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